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I | Abstract. In this paper we present new ESO/VLT FORS1 and ISAAC images of the fields around eight gravi- 

tationally lensed quasars: CTQ 414, HE 0230-2130, LBQS 1009-0252, B 1030+074, HE 1104-1805, B 1359+154, 
' H 1413+117 and HE 2149-2745. When available and deep enough, HST/WFPC2 data were also used to infer the 

photometric redshifts of the galaxies around the quasars. The search of galaxy overdensities in space and redshift, 

■ as well as a weak-shear analysis and a mass reconstruction are presented in this paper. We find that there are 
most probably galaxy groups towards CTQ 414, HE 0230-2130, B 1359+154, H 1413+117 and HE 2149-2745, 
with a mass <4 10 14 Mq1i -1 . Considering its photometric redshift, the galaxy group discovered in the field around 
HE 1104-1805 is associated with the quasar rather than with the lensing potential. 

Key words, galaxy clusters - gravitational lensing - individual HE 0230-2130, HE 1104-1805, B 1030+074, 
CL|' LBQS 1009-0252, B 1359+154, H 1413+117, CTQ 414, HE 2149-2745 
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, 1. Introduction ated with a galaxy cluster or a galaxy group. Actually, the 

■ . . . first X-ray detection of galaxy groups associated with the 
• • ■ Gravitational lensmg of distant quasars is a powerful tool , . , . , . n , , , , 

. i-ii i-i lensing potential m quadruple quasars has been made re- 

to address important cosmological and astrophysical ques- cently towardf , thg RXJ om+0551 (Morgan et aL 

tions. These include the distribution of dark matter m Q 095?+561 (chartas et aL 2002)) pQ m5+80 and 

£ . galaxies (Keeton, Kochanek & Falco 1998), the determi- B 1422+231 (Grant et al 2003) Optically, many galaxy 

nation of the Hubble parameter H through the mea- , , , . n . f . , , 

. *" , • clusters or groups have been discovered m the line of sight 

surement of the time delay between the quasar image , , ., ,. , , ,., _ . „„ 

,. ii i-ii towards gravitationally lensed quasars, like PG 1115+80 

light-curves (Refsdal 1964 , as well as ultimately the mea- , v , , lnn -N ^ mc - . CC1 , » . „,.„ . , , 

b r , i r i i-i a (Kundic et al. 1997), Q 0957+561 (Angonm-Willame et al. 

surement of the value of the cosmological constant A 19g4) thfl cioverleaf (Kneib et al. 1998b), RXJ 0911+0551 

(Kochanek 1996). All the complexity m resolving these (Kneib ^ &{ 2QQQ) ^ d gBS 1520+530 (Faure et al 2Q02) 

questions is due to the difficulty m modeling the lens- T ,, , mrrmnnn, , 

. • i tt -i r i i • i i ^ n this paper we present the new VLf/i'ORbl and 

me potential. Until now, most of the multiply imaged TC . . ^ , ,. , a u , . f . , , 

° * ' . .. . ISAAC observations and wide-field analysis of eight grav- 

quasars were thought to be lensed by one mam and single ,. , , , nmn A1 . 

, iiiii itationally lensed quasars m the southern sky: Gf Q 414, 

lensing galaxy, disregarding the local environment Yet, HE 0230 _ 2130; LBQS 1009 _ 0252; B i 03 0+074, HE 1104- 

while modeling these systems, it has been realized that 1805; B 135g+154; H 1413+n7 and HE 2149-2745. In pre- 

the addition of an external shear contribution is manda- . , ,. , r ,, n 

, vious studies of these systems, the modeling oi the lensing 
tory m many cases to reproduce the observations (Keeton, ;,-irrjr liri. t - i i. • 

,,„„... K , potential suffered from a lack of observational constraints 

Kochanek & Seliak 1997). It is now known, that m at least , , , n r 

J ' and a poor knowledge of tne quasar image surroundings, 

one fourth of the cases, the lensmg potential is much more r^, , ,. , , . , , , , . , , , 

' ° .... . I he modeling led indeed to the conclusion that the con- 
complex than a single lensing galaxy and that it is associ- , ., , . r , , , , r 

00 J tribution of an external shear was mandatory for most 

* Based on observations obtained with VLT/ANTU at ESO- of these cases to reproduce the quasar image configura- 

Paranal Observatory (programs 67.A-0502 and 69.A-0247) and tion using a realistic model for the main lensing galaxy. 

with the Hubble Space Telescope, operated by NASA. 
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In this paper, we first provide a brief summary of former 
studies for each system (Section |3J). Then we describe the 
new ground-based VLT/UT1 observations, and the HST 
dataset retrieved from the archive (Section EJ. In Section 
0]we discuss and give the results of our search for galaxy 
overdensities, and in Section [S] we provide the galaxy pho- 
tometric redshifts. In Sectional we present the weak-shear 
analyses and the mass reconstructions performed from the 
deep FORS1 images. And finally, in Section we pro- 
vide final remarks and present our conclusions. We adopt 
#o=65 km s- 1 Mpc" 1 , fi=0.3 and A=0.7 throughout the 
paper. 



2. Lens sample 

The goal of this project is to map the line of sight towards 
a large number of lensed quasars in the southern sky, up 
to the redshift of the lens. Indeed, the systems studied 
in this paper have not yet been modeled in a satisfactory 
manner due to the lack of constraints on the line of sight 
distribution towards the quasar on large scale. In particu- 
lar, little information is available in order to characterize 
the external shear due to possible intervening galaxy clus- 
ters or groups at high redshift along their line of sight. 
In this section, we present a brief summary of the previous 
studies performed on these systems. The quasar coordi- 
nates and redshifts, as well as some characteristics of the 
lensing potential are given in Table Q (from CfA- Arizona 
Space Telescope Lens Survey 1 , CASTLES; Muhoz et al. 
1998). 

As in most papers about gravitationally lensed quasars, 
we will use the convention that the brightest quasar im- 
age (at the date and in the filter of the discovery) is labeled 
A, the second bright one, is labeled B, etc. 

2.1. CTQ 414 

CTQ 414 (also named Q J0158-4325) has been first ob- 
served in the Calan-Tololo Survey (CTS; Maza et al. 
1995), and was fully identified as a doubly imaged quasar 
at redshift z =1.29 by Morgan et al. (1999). The quasar 
images are separated by 1.22". A PSF subtraction of the 
quasar images reveals an extended object lying in between 
the quasar images: it is thought to be the lensing galaxy. It 
is located 0.72" West to image A. The presence of the lens- 
ing galaxy has been confirmed by HST NICMOS obser- 
vations (CASTLES). No mass model has been performed 
for this system until now. 

2.2. HE 0230-2130 

HE 0230-2130 has been discovered in the course of the 
Hamburg/ESO survey (Wisotzki et al. 1996), and con- 
firmed to be a lensed quasar by Wisotzki et al. (1999). 
It is a quadruply imaged quasar at z =2.162, with image 



1 Gravitational Lens 
http:/ /cfa-www.harvard.edu /castles/| 



Database: 



angular separations from 0.7" to 2.1". The lensing galaxy 
was also observed and its apparent magnitude suggests 
that it is at redshift z; <1.6. No mass model has been 
proposed so far. 

2.3. LBQS 1009-0252 

LBQS 1009-0252 was studied by Surdej et al. (1993) and 
Hewett et al. (1994). It is a doubly imaged quasar at red- 
shift z =2.74. The images are separated by 1.53" and 
show strong Mg II absorption lines at z =0.869. A second 
quasar, quasar C, at redshift z =1.62, lies by chance at 
only 4.6" from quasar image A. The lensing galaxy was 
detected, close to quasar image B, in the HST I- and V- 
band frames (Lehar et al. 2000). 

In the lens model of Lehar et al. (2000), the host galaxy 
of quasar C is supposed to produce a shear of 7 ^0.053 
(after correcting for redshift differences) and, according 
to these authors, could dominate the external shear when 
fitting the main lensing galaxy by a Singular Isothermal 
Ellipsoid (SIE). A more recent model of the lens , requires 
an external shear: 7=0.017±0.009, with an orientation: 
PA 7 = 11±16 deg, if the lensing galaxy is modeled by a 
Singular Isothermal Sphere (SIS)(Claeskens et al. 2001). 



2.4. B 1030+074 

The bright radio source B 1030+074 is a two-component 
gravitationally lensed system at z =1.535 (Xanthopoulos 
et al. 1998). It was discovered during the course of the 
Jodrcll-Bank VLA Astrometric Survey (JVAS, Patnaik et 
al. 1992, Patnaik 1993, Browne et al. 1998, Wilkinson et 
al. 1998). The two quasar images are separated by 1.56". 
The lensing galaxy was discovered by Fassnacht & Cohen 
1998, at a redshift z;=0.599. It shows substructures that 
could be intrinsic to the galaxy or that could trace in fact 
an interacting galaxy system (Jackson et al. 2000). 
A first model of the lensing potential was proposed by 
Xanthopoulos et al. (1998). They modeled the lensing 
galaxy by a SIE and a flux ratio between the quasar com- 
ponents of 20:1. They have estimated the mass of the lens- 
ing galaxy within the Einstein ring to be 1.3 x 10 11 M©. 
Then, Lehar et al. (2000) proposed a two-component lens 
model, the main lensing galaxy G being modeled by a SIE, 
and the second one, C, modeled by a SIS. They concluded 
that, without external shear, the orientation of the SIE is 
inconsistent with that observed for the light of galaxy G, 
whereas it provides a good fit of the lens constraints. This 
disagreement is still pending. 

2.5. HE 1104-1805 

HE 1104-1805 (Wisotzki et al. 1993) is a doubly imaged 
quasar at redshift z =2.319, with image angular separa- 
tion of 3.19". In 1998, Courbin et al. discovered the lens- 
ing galaxy from ground-based near-IR (J- and X-band 
observations). Remy et al. (1998) have also detected the 
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Table 1. Properties of the lensed quasar sample. Column 1: quasar name. Columns 2, 3: right ascension and declination 
of quasar image A (J2000). Column 4: quasar redshift. Column 5: lensing galaxy redshift, where "?" means "unknown" 
and "(p)" means "photometric redshift" . Column 6: size of the system corresponding to the maximal angular separation 
between the quasar images (in arcsecond). Column 7: system configuration. Column 8: predicted value of the external 
shear strength and position angle of the shear (1): Lehar et al. 2002, (2): Saha & Williams 2003, (3): Lopez et al. 1998, 
(4) Rusin et al. 2001. The symbol "-" means that the system has never been modeled. 



Quasar 


RA 


DEC 


Zqso 




Size 


Configuration 


Ext. shear and 


PA(deg) 


CTQ 414 


01 ft 58 m 41.43 ! 


-43°25'3.4" 


1.29 


? 


1.22" 


double 




HE 0230-2130 


02 h 32 m 33.P 


-21°17'26" 


2.16 


<1.6 


2.10" 


quadruple 




LBQS 1009-0252 


10 fe 12 m 15.71 s 


-03°07'02" 


2.74 


0.88(p) 


1.54" 


double 


0.053 1 


-36 


B 1030+074 


10"33 m 34.08 ! 


+07°11'25.5" 


1.54 


0.599 


1.56" 


double 


0.09 1 


-39 


HE 1104-1805 


ll h 06 m 33.45 s 


-18°21'24.2" 


2.32 


0.73 


3.19" 


double 


0.05 1 


-38 


B 1359+154 


14 h 01 m 35.55 s 


+15°13'25.6" 


3.235 


? 


1.71" 


sextuple 


0.12 to 0.26 4 


71 to 97 


H 1413+117 


14 ft 15 m 46.40 s 


+ 11°29'41.4" 


2.55 


0.9(p) 


1.10" 


quadruple 


0.08^ 


70 


HE 2149-2745 


2i"52 m 07.44 s 


-27°31'50.2" 


2.03 


0.5(p) 


1.71" 


double 


>0.21 a 


58 



lensing galaxy in their WFPC2/7-band image, but not 
in their V-band image. Therefore, they suggested that 
the lensing galaxy is an early type galaxy. The redshift 
of the lensing galaxy is z;=0.729 (Lidman et al. 2000). 
The recent time-delay measurement (Ofek & Maoz 2003) 
gives a value of Ai=-161±7 days (68% confidence level), 
in disagreement with former results by Gil-Merino et al. 
2002 (Ai=-310±20 days) and Pelt et al. 2002 (Ai=-255 
days). Moreover, this system is known to show conspic- 
uous microlcnsing effects (Wisotzki et al. 1993, Courbin 
et al. 2000, Gil-Merino et al. 2002, Schcchtcr et al. 2003), 
making it more difficult to model. Several mass models 
have been performed (Wisotzki et al. 1998, Courbin et al. 
2000, Schechter et al. 2003), all leading to the conclusion 
that the contribution of an external shear is mandatory 
to reproduce simultaneously the quasar image configura- 
tion, the time-delay and the flux ratios, using a realistic 
cllipticity and mass for the lensing galaxy. 

2.6. B 1359+154 

The complex lensed system B 1359+154 consists in a set 
of six images of a quasar at z =3.235 (Rusin et al. 2001). 
The main lensing potential is due to a group of three galax- 
ies at z~l. Rusin et al. have assumed that the true number 
of lensed quasar image's is nine, and that only six of them 
are amplified and detected. They derived mass models, as- 
suming four different galaxy potentials for the three main 
lensing galaxies: SIS, SIE, Pscudo-Jaffc spherical (PJS) 
and Pscudo-Jaffc ellipticals (PJE) (for the definitions, sec 
Keeton & Kochanek 1998). None of these models matches 
well the flux ratios and the relative distances between the 
quasar images. The complexity of this system requires 
more observational constraints to be correctly modeled. 

2.7. H 1413+117 

The Cloverleaf is a distant quasar at z =2.558 discov- 
ered by Hazard et al. (1984) and identified as a lensed 



quasar by Magain et al. (1988). H 1413+117 shows four 
images, with angular separations from 0.76" to 1.1". The 
lensing galaxy was discovered by Kneib et al. (1998b) in 
the NICMOS-2/F160w observations. A galaxy overdensity 
was also discovered on the WFPC2/7- and i?-band (fil- 
ter I: F814w, filter R: F702w) images, in the magnitude 
ranges 23<I<25 and 23<R<25 (Kneib et al. 1998a,b). 
According to Kneib et al. (1998a), the color of the galax- 
ies building the overdensity suggest a redshift of z =0.9. 
Several models have been performed for the lensing object, 
using also constraints from the CO-map of the quasar im- 
ages (Kneib et al. 1998b), and the ones which reproduce 
best the observations are those taking into account the 
galaxy overdensity discovered by Kneib et al. (1998a, b). 

2.8. HE 2149-2745 

The doubly imaged quasar HE 2149-2745 at redshift z = 
2.033, presents an angular separation of 1.71" (Wisotzki 
et al. 1998). The lensing galaxy, appearing after PSF sub- 
traction of the quasar images, has an i?-band magnitude 
R=20.4. The V-R color of the lensing galaxy, if elliptical, 
corresponds to a maximum redshift of 0.5. The lensing 
galaxy lies at almost equal projected distance from both 
quasar images, and these images fall on a direction coin- 
cident with the major axis of the lensing galaxy. Lopez 
et al. (1998) proposed a model for the lens. They used a 
singular isothermal elliptical mass distribution (SIEMD, 
Kassiola & Kovner 1993). Due to the singular geometri- 
cal configuration of the quasar images, an external shear 
contribution is mandatory in the model. Because of the 
non-detection of any galaxy cluster in the field (mag- 
nitude limit: R^23), Lopez et al. imposed the external 
shear 7 to be less than 0.21. Their best models give a 
lens mass of : M=1.5(2.4) x 10 n M© for a lensing galaxy 
at redshift zj=0.3(0.5), corresponding to a mass to light 
ratio of M/£r=5(2). A time delay between the quasar 
images light-curves has been measured: A£=103±12 days 
(Burud et al. 2002). 
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Table 2. F0RS1 and ISAAC observations. Column 1: 
quasar name. Column 2: filter. Column 3: total exposure 
time in kilo-seconds. Column 4: field of view used for the 
analysis. 



Quasar 


Band 


Seeing 


Exp Time 


Field 


CTQ 414 


R 


0.50" 


1.45 


6.5' x6. 6' 




J 


0.51" 


2.00 


2.3'x2.4' 




Ks 


0.50" 


1.50 


2.2' x2. 2' 


HE 0230-2130 


R 


0.74" 


3.40 


6.4' x6. 6' 




J 


0.58" 


2.60 


2.5'x2.4' 




Ks 


0.50" 


2.70 


2.l'x2.4' 


LBQS 1009-0252 


R 


0.72" 


3.90 


6.5' x6. 6' 




J 


0.44" 


3.80 


2.4' x2. 5' 




K 


0.53" 


4.50 


2.2' x2. 2' 


B 1030+074 


R 


0.48" 


6.00 


6.6' x6. 6' 




J 


0.60" 


3.60 


2.4' x 2.4' 




Ks 


0.60" 


4.80 


2.3' x2. 3' 


HE 1104-1805 


R 


0.66" 


0.80 


6.6' x6. 7' 




J 


0.51" 


3.50 


2.4' x 2.4' 




Ks 


0.34" 


4.80 


2.2' x2. 3' 


B 1359+154 


R 


0.62" 


4.80 


6.6'x6.7' 




J 


0.64" 


3.60 


2.3' x2. 3' 




Ks 


0.47" 


4.50 


2.4' x 2.4' 


H 1413+117 


R 


0.58" 


1.70 


6.2' x6. 2' 




J 


0.59" 


2.60 


2.6' x2. 6' 




Ks 


0.49" 


3.60 


2.6' x2. 6' 


HE 2149-2745 


R 


0.57" 


1.70 


6.5'x6.7' 




J 


0.48" 


2.60 


2.3' x2. 5' 




K 


0.47" 


3.60 


2.3' x2. 5' 



3. Observations - data reduction 

3.1. Ground-based dataset 

We took deep and high resolution FORS1 i?-band images 
in order to exploit wide field images in the analysis 
of the lensing potential at large scale, performing a 
galaxy overdensity search and a weak-shear analysis. The 
near-IR images taken with the ISAAC instrument where 
mandatory to trace the stellar mass of galaxies at high 
redshift (z^3-4), and to measure galaxy photometric 
rcdshifts from multi-color analysis. 

The ground-based observations consist in two datasets 
taken with ISAAC (Cuby et al. 2002) and with 
FORS1 (Szeifert 2002) attached to the 8.2 m telescope 
VLT/ANTU, at ESO/Paranal. The imaging datasets were 
obtained during ESO periods 67 and 69 (programs 67. A- 
0502(A,B,C) and 69.A-0247(A,B,C)). The ISAAC images 
(1024 x 1024 pixels, 0.1484" per pixel) were taken in two 
near-IR filters: the J- and K s-filters. We reduced the data 
following standard procedures (ISAAC Data Reduction 
Guide, Amico et al. 2002). All the images (J- and Ks- 
bands) taken in P67 were affected by the "odd-even col- 
umn effect" , arising from an ISAAC mis-functioning dur- 
ing this period (this concerns objects CTQ 414, HE 0230- 
2130, H 1413+117 and HE 2149-2745). This artifact was 



Table 3. WFPC2 observations. Column 1: quasar name. 
Column 2: filter. Column 3: proposal P.I. and proposal ID. 
Column 4: total exposure time (in ks). 



Quasar 


Band 


P.I and I.D. 


Exp Time 


CTQ 414 


F555w 


Falco (9133) 


1 .60 




F814w 


Impey (8268) 


2.00 


HE 1104-1805 


F814w 


Surdej (5958) 


1.00 


B 1359+154 


F555w 


Impey (8268) 


5.20 




F814w 


Impey (8268) 


5.00 


H 1413+117 


F814w 


Turnshck (5621) 


5.20 


HE 2149-2745 


F555w 


Falco (8175) 


1.60 




F814w 


Falco (7495) 


1.60 



successfully removed through a Fourier analysis procedure 
(see ISAAC Handbook, Cuby et al. 2002). When images of 
a same object were taken during different nights we have 
reduced all images to the same "seeing", through convo- 
lution by a Gaussian, before combining them. A summary 
of the observations is given in Table [21 
The FORS1 i?-band images (2040 x 2040 pixels, 0.2" 
per pixel) were reduced using sky flat-fields and standard 
IRAF procedures for optical image reduction. Saturated 
stars were masked in order to perform the sky subtraction 
and the galaxy extraction (done both using SExtractor; 
Bcrtin & Arnouts 1996) . Quantitative information on the 
quality of the images is provided in Table [2] 
The sizes and positions of the ISAAC fields relative to the 
FORS1 fields are shown on FigOJand Fig§] 
As shown in Table the total exposure time varied from 
one system to the other, a consequence of the observations 
being conducted in service mode. In particular, some sys- 
tems are under-exposed like HE 1104-1805 in the i?-band 
and CTQ 414 in the Ks-bemd. 

3.2. HST imaging 

We have searched the HST public archive for extra im- 
ages of the quasar fields around: CTQ 414, HE 0230-2130, 
HE 1104-1805, B 1359+154, H 1413+117 and HE 2149- 
2745, in the 7- and ^-bands (F814w and F555w). Details 
of the observations are given in Table |3J In most cases we 
are able to use directly the WFPC2 association frames 2 , 
otherwise we used standard reduction procedures from the 
IRAF/STSDAS package. To project the HST images onto 
the sky we used the SWarp package from Bertin 3 . 

4. Looking for galaxy overdensities 

4.1. Photometry and galaxy-star separation 

We have searched for galaxy number overdensities in the 
field around the quasar images. Indeed, if these overden- 
sities are tracing galaxy clusters or galaxy groups, they 
must take part in the total lensing potential. 

2 http:/ /archive. eso.org/wdb/wdb/hst/asn_wfpc2/form 

3 http:/ /terapix. iap.fr/complement/oldSite/soft/swarp/ 
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3.5 p 




3.0 




3.0 




Fig. 1. Logarithm of the number of galaxies as a function of magnitude. Top panel: in the FORS1 (i?-band) fields 
around the lensed quasars and in the FORS Deep Field (FDF) Survey (thick solid line; Heidt et al. 2003). All the 
plots are rescaled to a 6.6' x 6.6' field of view. Middle panel: in the ISAAC (Ks-b&nd) fields around the lensed quasars. 
Bottom panel: in the ISAAC (J-band) fields around the lensed quasars. Diamonds: B 1359+154; squares: H 1413+117; 
long-dash line: HE 2149-2745; solid line: HE 1104-1805; dotted line: LBQS 1009-0252; dash line: B 1030+074; dash-dot 
line: CTQ 414; dash-dot-dot-dot line: HE 0230-2130. 
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Fig. 2. ISAAC ATs-band images of HE 0230-2130 (left) and H 1413+117 (right). The contours (from 2a, increasing by 
lcr steps) outline regions of same galaxy number density (label- value attached to each contour in galaxies per arcmin 2 ). 
North is to the top, East to the left. For HE 0230-2130, the galaxy overdensity is located Aa=40",A<5=-3" relative 
to quasar image A (coordinates (0,0)) and occurs in the magnitude range 19<Ks<22. For H 1413+117, the galaxy 
overdensity is located Aa=-3" /A6—-6" relative to quasar image A (coordinates (0,0)) in the same magnitude range 
as for HE 0230-2130. 



Table 4. Characteristics of the galaxy number overdensities. Column 1: quasar name. Column 2: limit for the galaxy 
catalog (said to be as complete as the FDF galaxy catalog where the curves in Fig^top do join). Column 3: magnitude 
range over which we can study the significance of the galaxy overdensity. Column 4: position of the center of the galaxy 
number overdensity relative to quasar image A (see Figs and QJ- The uncertainty on the exact position of the galaxy 
overdensity is about 5". Column 5: galaxy density at the location of the galaxy overdensity center. Column 6: FDF 
mean galaxy density in the magnitude range given in column 3. Column 7: FDF maximal galaxy density in the same 
magnitude range. Column 8: detection level (in a) of the cluster compared to the FDF mean galaxy density level. 
Column 9: detection level of the cluster (in the magnitude range given in parentheses) above the galaxy background 
level measured in the same frame. 





Compl. 


Mag- 


Galaxy overdensity 


Density 


FDF 


FDF 


Above 


Above 




limit 


range 


location 




mean 


max 


FDF 
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The extraction and photometry of the objects from 
all frames were performed using the SExtractor soft- 
ware (Bertin & Arnouts 1996). Galaxies are sepa- 
rated from stars using a combination of FWHM (Full 
Width-Half Maximum) versus magnitude (MAG_BEST in 
SExtractor) analysis and [imax (peak surface brightness in 
mag arcsec~ 2 ) versus magnitude analysis (see Bardeau et 
al. 2004). The stars are the brightest objects of the cata- 
logs, having the highest peak surface brightness and the 
smallest FWHM. 

4.2. Overdensities above the galaxy background level 

To search for galaxy number overdensities, we have 
proceeded as follow: for each field, the galaxy magnitude 
catalogs were split in magnitude slices of one magnitude 
width, allowing a first-order study of different redshift 
slices. For each magnitude slice, we have cut the image 
in a grid of 10x10 cells 4 . For each galaxy G\ in a cell 
we have measured the distance between G\ and its ninth 
closest galaxy Giq. Then, we have calculated the area, A, 
covered by a circle of radius 1Z=\Gi — G\q\. By definition, 
the number of galaxies included in the area A, with 

4 ~ 40" size cell in the /2-band image, 15" in the J- and 
Ks-im&ges 



radius TZ—\G\ — Giq\, is 10. Therefore, the galaxy number 
density associated to the galaxy G\ is the ratio 10/.4 
. A mean galaxy density for each cell is then derived, 
averaging the densities associated to all galaxies in the 
cell. The grid is chosen considering that the central part 
of a galaxy cluster would cover between 1 and 4 cells, 
depending on its redshift and on its central density. Then, 
the isodensity contours of the number of galaxies per 
arcmin 2 across the field are drawn: when an overdensity 
is observed in two consecutive magnitude slices, the two 
sub-catalogs are added in order to increase the signal 
to noise of the overdensity detection. The results are 
displayed in Fig. [3] and Fig.0] 



To evaluate the significance of these detections we 
have performed the following tests. As mentioned above, 
a galaxy cluster should appear in one to four contiguous 
over-densed cells in comparison to the rest of the frame. 
We have plotted the number of cells versus the galaxy 
density in the cells. This plot can be fitted by a Gaussian 
(see Faure et al. 2003 for a similar analysis), providing 
the mean galaxy density in a field for a given magnitude 
range, as well as a standard deviation a (Fig. Eland Fig. 
BJ. This estimator of detection is given in column 9 in 
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Tableland results are discussed in Section l4~4l 



4.3. Comparison with the FORS Deep Field 

For the i?-band images we have also used the FORS Deep 
Field 5 (FDF; Heidt et al. 2003) i?-band galaxy catalog as 
a reference field. The i?-band data set of the FDF covers 
an area of ~7.4'x7.4'(26.4 kilo-second exposure time), up 
to R=26.7 (50% completeness). The FDF field is centered 
close to the South Galactic Pole (l 7l 6 m 3.6 s , -25°45'46", 
J2000). It has been chosen because there was a priori no 
previously known galaxy cluster in this field. 

The FDF dataset is much deeper than our dataset 
(see Fig.nj. So we are allowed to compare them up to the 
magnitude corresponding to the 100% completeness limit 
of our dataset. These magnitudes are displayed in Tabled 

We have proceeded as follows: first, we have split the 
FDF catalog in magnitude ranges identical to those where 
galaxy ovcrdcnsitics were discovered in the fields around 
the quasars. When an overdensity was detected in a mag- 
nitude range such that the upper magnitude was higher 
than the completeness magnitude, we have compared only 
the part of the catalog for which the upper magnitude 
was equal to the completeness magnitude (see Table 21 
columns 2 and 3). We have partitioned the FDF in a grid 
of 10 x 10 cells. For each cell, a galaxy density in a given 
magnitude range is measured and as detailed in Section 
14.21 we have fitted a Gaussian to the histogram describing 
the "number of cells" versus the "galaxy density" . This 
provides the mean galaxy density, the standard deviation 
and the maximal density expected in a cluster-free region, 
in each magnitude range. The results of the comparison 
between the FDF and our dataset are displayed in Tabled 



4.4. Results of the search for galaxy overdensities 

The galaxy overdensities detected towards CTQ 414, 
B 1030+074, HE 1104-1805 and H 1413+117 are well 
above the FDF mean level, revealing a 2 dimensional 
structure in the line of sight towards these quasars. 
On the contrary, the galaxy overdensities discovered 
towards HE 0230-2130, LBQS 1009-0252, B 1359+154 
and HE 2149-2745 show values lower than (or very close 
to) the maximal density measured in the FDF field. 
This could mean that these densities are compatible 
with the variance of the counts within the field of view. 
Hence, they may correspond to areas where the count 
fluctuations reach a slightly higher level, still compatible 
with random deviation from the mean density. Another 
possibility is that there are undetected galaxy groups 
in the FDF, leading to values of the galaxy densities as 

5 The FDF photometric catalog is available from: 
http: / /cdsweb.u-strasbg.fr/cgi-bin /qcat? J / A+ A/398 /49 



strong as in the lens fields. 

Comparing the detected overdensity to the mean galaxy 
density measured over the rest of the frame, we see that, 
apart from the galaxy overdensities towards HE 0230- 
2130 and LBQS 1009-0252 (respectively detected at 2.5ct 
and 3.5cr above the background level), the overdensities 
measured towards the other six lensed quasars are all 
significant (4<r to 6er). Notice that the high detection level 
(26<t) of the cluster towards HE 1104-1805 compared 
to the FDF galaxy number density is most probably 
overestimated. This results from the fact that the galaxy 
densities are analyzed over a very narrow magnitude 
range (Am=0.3), and hence, a small galaxy sample 
(^100 galaxies), leading consequently to possibly large 
statistical effects. 

Furthermore, there are two cases for which we have 
detected more than one cluster in the line of sight 
towards the quasar images. Towards H 1413+117 a 
second galaxy overdensity has been detected in the 
magnitude range 25<R<28, at 3er above the background 
level (centered 36" South- West to the quasar images). 
In a similar way, we find a second overdensity towards 
HE 1104-1805 (78"South-West) in the magnitude range 
24<R<25, at 5er above the background level. But in these 
two cases, the data used for the estimation extend below 
the catalog completeness limit, and therefore we cannot 
assure whether these are real detections or whether they 
correspond to standard deviations at that level. 

In most cases, the number of galaxies detected in the 
near-IR frames is too low to perform a secure iso-density 
analysis. Nevertheless, in the Ks-band images towards 
HE 0230-2130 and H 1413+117 we find again the galaxy 
overdensities seen in the i?-band images (see Fig. [3J| . The 
shifts between the positions of the overdensities in the 
R- and Ks- band images come mainly from the use of 
different cell sizes and locations. 

According to these results we conclude that, apart 
from the marginal galaxy overdensity towards LBQS 1009- 
0252, all others are possible 3D-structures. That is what 
we are going to check in the following sections. 

5. Photometric redshifts and color-magnitude 
analysis 

We can test whether the galaxy overdensities discovered 
in Section 0] are associated to galaxy clusters or galaxy 
groups by measuring the galaxy photometric redshifts. 
The redshift will also tell us whether these structures lie 
in the line of sight towards the quasar or whether they are 
at the quasar redshift. 

5.1. Photometric redshifts 

Whenever a galaxy overdensity was discovered on the 
line of sight towards the quasar images, we performed the 
photometric redshift analysis using the Hyperz software 
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(Bolzonella et al. 2000) and the full dataset available 
(described in Section The input to Hyperz are 
the galaxy magnitudes in at least 4 filters 6 . Hyperz 
performs a comparison between the photometric spectral 
energy distribution (SED) of the observed galaxy and 
those derived from a set of reference template spectra, 
using the same photometric system. Reddening is taken 
into account using the Calzetti law (Calzetti et al. 2000). 
We estimate that the current dataset is deep enough to 
reach giant elliptical galaxies at redshift up to z =3.5. 
Two parameters are important to test the significance 
of the photometric redshifts found: the x 2 value of the 
SED fit, and the probability of obtaining this x 2 : P(x 2 )- 
In the following analysis, we only take into account 
the redshifts obtained with x 2 ~1 and P(x 2 )>70%, 
considering that this is the limit of confidence for the 
redshift determination. 

In the FORS1 field around CTQ 414 we have 
computed the galaxy photometric redshifts using the 
VLT and HST data. The overlapping field covers only a 
very small area (see Fig. [SJ, hence only 7 galaxies have 
their redshift determined using 5 magnitudes (V, R, I, 
J and Ks). In addition, we tentatively determined the 
redshifts of 12 galaxies from their V-, R-, J- and i^s-band 
magnitudes, and of 9 galaxies from their R-, I-, J- and 
Ks-band magnitudes. The results are plotted in Fig. [SJ 
showing a significant peak in redshift at z =0.5 ±0.1. In 
the following, this is assumed to be the redshift of the 
galaxy overdensity discovered towards the quasar. 

For HE 0230-2130, we have not been able to deter- 
mine the photometric redshifts of the galaxies building 
the observed marginal overdensity, because of the low 
depth of the WFPC2 images available for this field. 
However, according to Wisotzki et al. (1999) the lensing 
galaxy redshift is z; <1.6. Assuming that the lensing 
galaxy is part of the galaxy cluster tentatively discovered, 
we will use arbitrarily z =1.0 for the lensing potential 
redshift in the mass reconstruction. 

For LBQS 1009-0252 there is no significant galaxy 
overdensity in the field (Section and with the galaxy 
photometry available in only 3 bands (R, J and Ks) 
we cannot perform a photometric redshift analysis. 
Therefore, a weak-lensing analysis is necessary to define 
whether or not there is a dark source of shear generating 
the external-shear value predicted by Lehar et al. (2000, 
see Table [TJ. 

For B 1030+074, the lack of multi-wavelength data 
makes impossible photometric redshift measurements. 
Moreover, the galaxy overdensity detected in the FORS1 
field around B 1030+074 (see Section 0} is centered 

6 A study of the accuracy of the code, according to the num- 
ber of filters used to infer the redshifts, is detailed in Bolzonella 
et al. (2000) 



outside the ISAAC field. Would this galaxy overdensity 
be a cluster or a group, it would be too far away from the 
quasar images (~2') to play any significant role in the 
lensing potential (shear strength negligible). Hence, as 
for LBQS 1009-0252, we expect the weak-shear analysis 
to conclude whether or not there is a galaxy group or 
cluster in the line of sight towards the quasar. 

Around HE 1104-1805 we have derived the photo- 
metric redshifts of the galaxies present in the overlapping 
HST and VLT fields (see Fig. 0J. One WFPC2 chip 
covers mainly the galaxy overdensity located 50" North 
to the quasar images. The four galaxies detected both on 
this chip and on the three VLT images have a redshift in 
the range: 1.95<z<2.10. Therefore, this cluster is most 
probably associated with the quasar (z=2.33) rather than 
with the lens. The photometric redshift of the second 
galaxy overdensity has not been measured, being mostly 
outside the overlapping zone of the VLT and HST fields. 

In the field around the six-image lensed quasar 
B 1359+154 we have measured the photometric red- 
shifts of galaxies appearing in the overlapping VLT and 
HST fields (see Fig|SJ|. There are very few galaxies for 
which the SED is fitted with a x 2 ~1 (16 galaxies among 
the 34 detected in the four bands). Half of them have a 
redshift between 1.1 and 1.5 and might correspond to the 
galaxy overdensity detected towards B 1359+154 because 
of their location in the field. Therefore, we will consider 
in the following, that the redshift of the galaxy cluster 
towards B 1359+154 is z =1.3. 

In the field around H 1413+117 we have computed 
the photometric redshifts for galaxies located in the area 
covered by both the VLT and HST datasets. Two peaks 
are observed in redshift: a first one at z =0.8±0.3, and 
a second one, at z =1.75+0.2. As the two overdensities 
are mostly superimposed on the common HST and VLT 
frame, it is hard to identify their corresponding redshift. 
Yet, in the original paper reporting on the discovery of the 
South-East overdensity discovery (Kneib et al. 1998b), 
the photometric redshift derived (from three filters only) 
was z— 0.9+0.1. We will assume in the following that 
the redshift of the South-East overdensity is z=0.8±0.3. 
The second galaxy overdensity could correspond either to 
the absorption system at z=1.661 or to that at 2=1.87 
(Turnshek et al. 1998; Magain et al. 1998.) 

Finally, towards HE 2149-2745 the distribution of 
photometric redshifts (using the HST and VLT images) 
shows a peak in redshift at z =0.7+0.1 with a high 
confidence level. This suggests the presence of a galaxy 
cluster or galaxy group towards the doubly imaged 
quasar. On the other hand, previous studies by Wisotzki 
et al. (1998) and Lopez et al. (1998) argued that the 
lensing galaxy maximal redshift is z =0.5, according to 
its V-R color. Due to the uncertainty on the photometric 
redshift determination and on the redshift limit given by 
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Wisotzki et al., it is still uncertain whether the lensing 
galaxy is part of the galaxy cluster/group discovered in 
this study. 



5.2. Color versus magnitude diagrams 

We have analysed the R-Ks color versus Ks magnitude di- 
agrams of the galaxies detected in the overlapping FORS1 
and ISAAC fields. Unfortunately, no elliptical red se- 
quence appears in any of the field. The presence of such a 
sequence would argue in favor of the presence of a massive 
galaxy cluster towards the quasar. There are possibly two 
reasons for a non-detection: I) there are too few galaxies 
in the common R- and Ks-band catalogs (between 50 and 
200). Hence, the elliptical red sequence of galaxy clusters 
does not show up in a conspicuous manner. 2) the galaxy 
clusters or groups are not massive enough to have an el- 
liptical red sequence. 

6. Weak-shear analysis and mass reconstruction 

For each field, we have searched for the weak-shear sig- 
nature that the galaxy groups or clusters should imprint 
on background galaxies. According to the gravitational 
lensing theory, the presence of a galaxy cluster or a galaxy 
group in the field generates distortions and amplifications 
of the background galaxies. In practice, we measure the 
"statistical" distortions: the mean ellipticity and mean 
position angle of the faintest galaxies in cells of a given 
size. 



6.1. PSF correction 

We have performed the weak-shear analysis on the 
FORSl/i?-band images of the systems. First of all, we 
have to correct the galaxy ellipticities for the instrumental 
distortions across the field. Each night, the PSF variations 
across the field are measured from a sample of stars, 
using the Im2shape software (Bridle et al. 2002, Kneib 
et al. 2003). The ellipticity of the PSF varies typically 
from e=(a-b)/(a+b)=0.002 to 0.05 across the FORSI 
field. Then, the background galaxies are deconvolved in 
an analytical way, modeling them by 2D-functions and 
using the closest PSF. 



6.2. Distorsion measurement 

Once the galaxy images have been deconvolved and using 
only galaxies with ellipticities known with an accuracy 
better than 20%, we can compute mean galaxy ellipticities 
and mean galaxy orientations across the field. The field is 
divided into a grid of cells. In order to optimize the coher- 
ence of the shear-map pattern, we attribute to any galaxy 
within the cell an ellipticity and an orientation derived 
through averaging the ellipticities and the orientations of 



all galaxies within a radius r slightly larger than half the 
cell size. Then, within each cell, the attributed ellipticities 
and attributed orientations are averaged for all the galax- 
ies. These final "mean ellipticity" and "mean orientation" 
are allocated to the cell center. In a field with a massive 
cluster (CI 0024+1654, Hudelot et al. 2004, private 
communication) we have tested that, for the total number 
of galaxies used to perform the weak-shear analysis (be- 
tween 1000 and 2500 galaxies per field), a combination of 
a cell -40"x40", and a radius r of 0.01° to 0.02°, allows 
to extract the statistical information in a satisfactory way. 

None of the shear map of these systems reveals a signal 
strong enough to be spotted at first glance. However the 
shear map can be used to derive the lensing potential mass 
distribution as well as an upper limit of its mass. This is 
detailed in the following subsection. 

6.3. Mass reconstruction 

The galaxy mean ellipticities and mean orientations are 
given as input to the LensEnt2 software (Bridle et al. 1998, 
Marshall et al. 2002), a code which allows to built the 
mass distribution of a lens using the shear-map. The red- 
shift of the lens plane assumed for this reconstruction is 
that provided in Tabled (column 4). The source plane red- 
shift is derived from the photometric redshift of the FDF 
galaxies (E. Labbe, private communication, see Heidt et 
al. 2003 for the photometry): in a given magnitude range 
(displayed in Table EJ), we have averaged the FDF galaxy 
photometric redshifts and assigned this mean redshift to 
the galaxy sample. For all magnitude ranges considered 
here, the mean galaxy redshift is z ~ 1.4. The reconstruc- 
tion is performed in a frame 4 times larger than the frame 
where the weak-shear was analysed since the shear pat- 
tern eventually observed in the FORSI field could be due 
to the presence of a mass potential located outside the 
field of view. LensEnt2 uses a maximum entropy method 
to reconstruct the mass potential. The program simulates 
a synthetic population of unlensed background galaxies at 
a same redshift and in a field of the same size as the ob- 
served one (Bartelmann et al. 1996). Then LensEnt2 adds 
in the field a mass clump, with a compactness character- 
ized by its Intrinsic Correlation Function (ICF, Marshall 
et al. 2002), at the lens redshift, and analyses the similar- 
ity between this simulated shear pattern and the observed 
shear pattern in each field. 

The ICF of the mass clump is assumed to be Gaussian. 
For each system, we have performed a set of mass re- 
constructions using different ICF widths. The "evidence" 
(see Marshal et al. 2002) determined for each mass recon- 
struction provides a good discriminator of the ICF width, 
weighting the reality of the obtained mass reconstruction. 
The width of the ICF finally selected is given in Tabic [5j 
Typically, the mass bump which reveals a galaxy cluster 
in a reconstruction map corresponds to an ICF width be- 
tween 100" and 300". 
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Table 5. Results of the weak-shear analysis and of the mass reconstruction. Column 1: quasar name. Column 2: 
number of galaxies per arcmin 2 corresponding to the galaxy overdensity. Column 3: position of the galaxy overdensity 
center relative to the quasar images. Column 4: photometric redshift of the galaxy overdensity and corresponding error 
bar ("*": value inferred from previous paper, see text for references). Column 5: magnitude range for the background 
galaxies used for the weak-shear analysis . Column 6: number of galaxies used for the weak-shear analysis and the 
mass reconstruction. Column 7: ICF width of the mass reconstruction in arcsecond. Column 8: upper limit of the mass 
integrated from the mass reconstruction in a radius of 500 kpc and centered at the galaxy group position. Column 
9: galaxy cluster velocity dispersion corresponding to the mass given in Column 8. Column 10: shear orientation in 
degree. 
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Since the noise properties of maximum-entropy inversion 
methods are not easy to interpret, we have performed, 
for each system, 80 reconstructions starting with modified 
catalogs of galaxies (positions are retained but ellipticity 
directions are assigned at random). For each of the 80 
mocked catalogs, mass reconstruction is performed, and 
a set of mass bumps is obtained. Integrating over the 80 
catalogs, the mass bump distribution shows a number of 
peaks. Their median gives the noise level of this procedure. 
Results of the mass reconstructions towards the systems 
are displayed in Fig. H3 The contours outline regions of 
same signal to noise ratio. The signal to noise of the mass 
reconstruction is low at the location of the galaxy over- 
density (S/N<3 in all frames), and there are regions in 
the field with higher mass signal than at the position of 
the galaxy overdensity. The low signal at the overdensity 
location could have different origins: 1) there is actually 
no mass clump at this position 2) the mass inferred in the 
fields are too low to be detected by the reconstruction, 
and therefore the masses measured here are upper limit 
3) the largest structures in the field (real or artifacts due 
to the low number of background galaxies in some cases) 
are adding noise, and the signal level at the overdensity 
position is biased. Therefore, we investigate the limit of 
the LensEnt2 software performing the following test: for 
each system, we have lensed the background galaxies with 
a Pseudo Isothermal Elliptical Sphere (Kassiola & Kovner 
1993) of small ellipticity (e=(a 2 -b 2 )/(a 2 +b 2 )=0.05) and 
centered on the quasar image A. We have looked for the 
smallest mass potential that the software can reconstruct. 
This value varies in each field as the number of background 
galaxies is different. The lower mass values (in a radius of 
500kpc) are displayed in Table El (column3), as well as 



Table 6. Galaxy group masses. Column 1: quasar name. 
Column 2: mass integrated in the original mass map, in 
a radius of 500 kpc centered on the galaxy overdensity. 
Column 3: lower mass value reconstructible by LensEnt2. 
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the mass integrated in a same radius, but centered on the 
overdensity in the original mass map. 

6.4. Results 

The sizes of the mass clumps computed by the program, 
as well as the masses inferred, correspond to small galaxy 
clusters, or even galaxy groups. 

In Table |SJ the mass displayed are integrated over a 
radius of 500 kpc (typical size for a galaxy cluster) and 
centered on the galaxy overdensity. The velocity disper- 
sion is derived from the mass value, assuming that the 
mass distribution follows an Isothermal Sphere potential 
with a core radius of 30 kpc. 

We have chosen to rebuild the mass towards 
LBQS 1009-0252 even though no evidence for the 
presence of any galaxy cluster towards the quasar has 
been found. We choose arbitrarily z=l for the lensing 
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Fig. 7. Shear strength induced by the galaxy groups to- 
wards the different systems as a function of Ho . 



potential redshift and a galaxy magnitude range of 
R=[23,28[ for the measurement of distortions. The lack 
of mass concentration confirms that we do not detect 
any extra gravitational potential in the field around 
LBQS 1009-0252. 

From the results of the mass limit test displayed in 
Table we see that the only systems for which the mass 
inferred M have a physical relevancy are HE 2149-2745 
and HE 0230-2130, as they appear to be larger than the 
detection limit M[ ensei i. In the other cases the masses 
have to be taken as upper mass limits. In the case of 
LBQS 1009-0252 the mass measured is well below the 
detection limit, and this suggest that the overdensity 
detected here is not related to a mass potential. 

Finally, the LensTool software provides a map of the 
shear induced by a mass distribution according to the dis- 
tance to the potential center. Therefore, we can easily ex- 
tract the shear strengths that the galaxy groups induce at 
the quasar positions. These values are displayed in Fig. d 
versus the value of the Hubble constant. 

The shears computed here are derived from the masses 
of the galaxy groups. They are, as a consequence, shear up- 
per limit values in the same case of CTQ 414, B 1030+074, 
B 1359+154 and H 1413+117. We notice that they are 
in general lower than those predicted in former studies 
(displayed in Table ^J. Moreover the orientations are also 
different from the predicted ones. These discrepancies be- 
tween our results and previous ones might come from the 



fact that the predicted shear values and orientations had 
been formerly deduced by modeling the lensing galaxy as 
a sphere, that is to say without taking into account its ori- 
entation and ellipticity, while these last parameters dom- 
inate the shear strength and orientation in most cases. 

7. Conclusions 

We have studied the surroundings of the quasar images 
in eight gravitationally lensed systems in order to search 
for the presence of galaxy clusters or galaxy groups 
which could contribute to the "external shear" requested 
from previous studies. We have used three different and 
independent methods: the search of galaxy number over- 
densities, the search of a peak in the redshift distribution 
of galaxies -using photometric redshifts-, and finally a 
weak-shear analysis combined with mass reconstruction 
of the background galaxies. 

There is convincing evidence for the presence of a 
galaxy group in the line of sight towards the following 5 
lensed quasars: CTQ 414, HE 0230-2130, B 1359+154, 
H 1413+117 and HE 2149-2745. Towards LBQS 1009- 
0252 and B 1030+074 we have not been able to compute 
the photometric redshift of the overdensity, and, therefore, 
we cannot conclude about their 3D reality. The galaxy 
overdensity discovered in the field around HE 1104-1805 
seems to be associated with the quasar rather than with 
the lensing potential. 

The mass reconstruction leads to small mass clump 
sizes, related to galaxy groups. Moreover, the galaxy group 
masses and the shear they induce at the quasar position 
have to be regarded as upper limits, as we are using the 
reconstruction mass method at the edge of its applicabil- 
ity. 

The next step of this work is to incorporate the infor- 
mation derived from the wide-field analyses in the mod- 
eling of the lensing potential towards these eight lensed 
quasars. This will ultimately lead to an improved value for 
the Hubble constant, as derived from available time-delay 
measurements between the quasar image light-curves. 
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Fig. 3. FORS1 i?-band images of H 1413+117, B 1359+154, B 1030+074 and CTQ 414. North is to the top, East 
to the left. The contours (from 2cr, increasing by la) outline regions of same galaxy number density (label-value 
attached to each contour in galaxies per arcmin 2 ). The dash-line squares represent the position and relative size of 
the ISAAC field, while the solid-line squares represent the position, size and orientation of the WFPC2 field (in the 
case of CTQ 414 , the dash-point squares represent the position, size and orientation of the WFPC2 y-band field, 
different from the one of the WFPC2 /-band ones). For each system, quasar image A is located at coordinates (0,0). 
For H 1413+117, 2 galaxies overdensities are discovered: at Aa=-10",A<5 =-10" in the magnitude range 23<R<25, 
solid-line contours and at Aa=35",A<5 =-5" in the magnitude range 25<R<28, dash-line contours. For B 1359+154, 
a galaxy overdensity is found at Aa=-8",A<5 =3" (22<R<23). For B 1030+074, also one galaxy overdensity is 
detected (Aa=-90",A<5 =-70", 21<R<23). Finally, for CTQ 414, a galaxy overdensity is discovered in the magnitude 
24<R<25 at Aa=-10",A<5 =-5". 



Fig. 4. FORS1 R-band images of HE 0230-2130, HE 1104-1805, LBQS 1009-0252 and HE 2149-2745. The contours 
(from 2cr, increasing by la) outline regions of same galaxy number density (label-value attached to each contour in 
galaxies per arcmin 2 ). The dash- line squares represent the position and relative size of the ISAAC field, while the 
solid-line squares represent the position, size and orientation of the WFPC2 field. North is to the top, East to the 
left. For HE 0230-2130, one galaxy overdensity is discovered at Aa=40",A<5 =-5" (24<R<25). For HE 1104-1805, 
two galaxy overdensities are discovered: at Aa=0",A<5 =50", 23<R<24, solid-line contours and at Aa=60",A(5 =- 
50", 24<R<25, dash-line contours. For LBQS 1009-0252, a galaxy overdensity is found at Aa=60",A(5 =40" to 
quasar image A (23<R<25). For HE 2149-2745, the closest galaxy overdensity to the quasar images is centered at 
Aa=8",A<5 =-7"(23<R<25). 
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Fig. 5. Number of galaxies versus photometric rcdshifts for CTQ 414 (top left hand panel), for HE 1104-1805 
(top central panel), for B 1359+154 (top right hand panel), for H 1413+117 (bottom left hand panel) and for 
HE 2149-2745 (bottom right hand panel). For CTQ 414, the solid-line corresponds to the photometric redshift 
distribution determined from the: V, R, I, J and Ks galaxy magnitudes; the dotted-line corresponds to the redshift 
determination using the: R, I, J and Ks galaxy magnitudes; the dash-line corresponds to the redshift determination 
using the: V, R, J and Ks galaxy magnitudes. 



Fig. 6. Mass reconstruction superimposed on the FORS1 images. The contours outline regions of same signal-to-noise 
(dash-lines: S/N<1, solid-lines: S/N>1). From the top left hand corner to the bottom right hand corner: (a) CTQ 414, 
(b) HE 0230-2130, (c) LBQS 1009-0252, (d) B 1359+154 and (c) H 1413+117. 
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